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Spectral Functions of Quarks in Quark Matter^ 
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Abstract. We present a simple albeit self-consistent approach to the spectral function of light quarks in 
infinite quark matter. Relations between correlation functions and collision rates are used to calculate the 
spectral function in an iterative procedure. The quark interactions are described by the SU(2) Nambu- 
Jona-Lasinio model. Calculations were performed in the chirally restored phase at zero temperature. 
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1 Introduction 

It is well known that short-range correlations have influ- 
ence on the properties of nuclear matter and finite nu- 
clei; a substantial amount of high-momentum processes 
is contained in the nucleon spectral function. There have 
been many theoretical approaches trying to understand 
the short-range correlations. Of particular interest are the 
self-consistent calculations for the spectral function of nu- 
cleons in nuclear matter from Lehr et al. Q. The results 
of their simple model are in good agreement with sophis- 
ticated calculations from many-body theory. It is striking 
that their model is very successful in describing the influ- 
ence of short-range correlations on the properties of nu- 
clear matter using a simple pointlike nucleon interaction 
with a constant scattering amplitude. 

Motivated by the success of this model we have taken 
up the concept to investigate the properties of light quarks 
in infinite quark matter . It is our basic assumption that 
the properties of the spectral function are dominated by 
phase space effects and the overall strength of the interac- 
tion just like in the case of the nucleons. The detailed 
structure of the interaction should be relatively unim- 
portant as long as the relative symmetries are respected. 
We use relations between the spectral function and the 
collisional self-energies to construct a simple albeit self- 
consistent model. The spectral function can then be cal- 
culated in an iterative process beyond the quasiparticle 
approximation. The quark interactions are described by 
the Nambu- Jona-Lasinio (NJL) model ||] . It has the same 
symmetries as QCD and describes an effective pointlike 
interaction. 

As a first step the model has been applied to the sim- 
plest system, namely the chirally restored phase at zero 
temperature. Mean field effects were neglected, diquark 
condensates that arise in the color superconducting phase 



were not included. For the collisional self-energies only 
the lowest order contributions, the Born diagrams, were 
considered. In the numerical calculations the influence of 
the coupling strength and the chemical potential on the 
properties of the spectral function was investigated. 



2 The Model 

In this section we will briefly review our model. For more 
details we refer to §. The underlying Gr een's function 
formalism is discussed in much detail in |^,^. Note that 
current quark masses are neglected throughout this work 
and only systems in thermal equilibrium are considered. 
We use the Nambu- Jona-Lasinio model to describe the 
quark interactions in our approach. It is an effective inter- 
action model that was designed to resemble the symme- 
tries of QCD. The SU(2) Lagrangian is given by 



(1) 



where G is the coupling strength, independent of energy 
and momentum, and the are isospin Pauli matrices. 
Due to the pointlike interaction this model is nonrenor- 
malizable and a three momentum cutoff A is introduced. 
Currently we do not consider any extensions to the La- 
grangian that lead to color superconductivity, cf. ||^ . 
The correlation functions 

*.g>(l,l') = WlMl')), 
-zff<(l,l') = Wi'Mi)) 

are the fundamental elements of our model. In thermal 
equilibrium they are related to the spectral function A 
via the thermal Fermi distribution np{pQ): 
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Fig. 1. Examples for diagrams that are resummed in our cal- 
culation. The lines correspond to free propagators. 



Note that the Green's functions and the spectral function 
are matrices in spinor space. The single particle self-energy 
can be decomposed into a mean- field part Z""* and the 
collisional self-energies . Thus the retarded self-energy 
^rct given by: 

+ 0{t,-ty)[S>{l,l')-S<{\X)\. (4) 

The time local mean-field self-energy corresponds to 
the motion of particles in a potential. It is responsible for 
the dynamical generation of the constituent quark masses. 
The collisional self-energies contain the effects of particle 
decays and collisions in the medium. In lowest order they 
are given by the Born diagrams (cf. left diagram in Fig. 
0). One finds 



(5) 



where and g^ are full propagators and the integrals run 
over P2,P2, and p^. Eqs. (^) look like total collision rates 
and can be used to determine the width of the spectral 
function: 



r{p) = -2imi:-*(p) = i[i:>{p) ~ s<{p)]. 



(6) 



The real part of 17''°' is related to by a dispersion re- 
lation. If the width is known over the full energy range 
Rei7'''^' can be calculated dispersively. Using F and Rei7''°' 
the spectral function can be explicitly written as: 



A{p) = -2Im.g''°'(p) = -2Im- 



1 



(7) 



Eqs. (||),(||), and (||)-(^ form a set of equations describing 
a self-consistency problem. A direct solution is not easily 
possible. It is possible, however, to solve the problem it- 
eratively by starting from an initial guess for one of the 
quantities. In this way we can find a solution for the spec- 
tral function. The effect of self-consistency can be nicely 
illustrated in the language of Feynman diagrams. Using 
full propagators g^ that depend themselves on in eqs. 

means that we sum nonperturbatively over a whole 
class of diagrams. Some examples are shown in Fig. ^ 

Finally, we have to discuss the matrix structure of the 
spectral function in spinor space. To find the most general 
form A can be decomposed in terms of the 16 independent 
products of the 7 matrices. Demanding invariance under 
parity and time reversal symmetry we find for the rest 
frame of the medium (in thermal equilibrium) |^: 

^(P) =Ps(po,P^)+Po(po,P^)7° + Pv(po,P^)p-7, (8) 
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Fig. 2. The width Fo and the spectral function po of quarks 
at different momenta. Solid lines correspond to the usual NJL 
coupling strength, dashed lines to a coupling twice as large. 
Dotted lines have been obtained with a coupling four times 
larger than the usual value. 

where p is a unit vector in the momentum direction. Note 
that the density of states is given by po alone and that ps 
has to be zero in the chirally restored phase. It follows from 
eqs. (|^)-(0) that F and must have the same structure 
as p. 

Currently we apply two simplifications to the model. 
First, we restore chiral symmetry 'manually' by setting 
E"^^ to zero. This can be done for any density and tem- 
perature since m = is always a solution of the equation 
for the constituent quark mass However, one has to 
be aware that this might not be the thermodynamically 
favored phase when also a finite solution for m exists. 
Second, we neglect the real part of Z""^"' due to techni- 
cal reasons For nuclear matter this was a reasonable 
approximation To make sure that the effects of this vi- 
olation of analyticity are not significant we have checked 
that ReZ"'°' is small compared to cf. eq.(m). 



3 Results 

All our calculations were performed at zero temperature 
and in the chirally restored phase. The calculations were 
initialized with a constant width, Fq = 1 MeV and F^ = 0. 
Self-consistency was achieved after two iterations. First, 
we chose the quark matter density such that it is compara- 
ble to regular nuclear matter, pq,„ — S-pnm — 3-0.17fm~^. 
This yields a Fermi energy oi cop — 0.268 GeV. The cutoff 
A and the coupling constant G of the NJL model were 
chosen so that the model reproduces the known values |^ 
for the quark condensate and the pion coupling constant 

in vacuum. To investigate the influence of the coupling 
on the spectral function we did also calculations with two 
times and four times larger coupling strengths. In reality 
quark matter with a chemical potential oiuip — 0.268 GeV 
would not be in the chirally restored phase. Therefore 
we made additional calculations with a higher density of 
Pqm = 1.53 fm^"^, corresponding to ujp = 0.387 GeV. This 
chemical potential is well beyond the chiral phase transi- 
tion in the NJL model [|). 

Fig. I shows our results for the spectral function and 
its width at several momenta and coupling strengths us- 
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Fig. 3. The width /o and the spectral function po of quarks at 
different momenta. Solid lines correspond to ujf ~ 0.268 GeV, 
dashed lines to lof — 0.387 GeV. 
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Fig. 5. Series of diagrams corresponding to dynamically gen- 
erated mesons. 

as large the short-range correlations increase but still the 
depletion effect is below 1%. Only for the largest cou- 
plings we find a high momentum tail of a few percent, 
comparable to the case of nucleons. In Fig. |4|(b) it can be 
seen again that the effect of the larger chemical poten- 
tial is similar to an increased coupling. The depletion for 
ujf ~ 0.387 GeV grows by almost one order of magnitude 
compared to ojp = 0.268 GeV. 




p [GeV] 



Fig. 4. (a) Quark momentum distribution in quark matter 
for UJF = 0.268 GeV and three different couplings (see Fig. ^ for 
details), (b) Momentum distribution for two different chemical 
potentials (see Fig. ^ for details). 



ing UJF — 0.268 GeV. Due to the pointlike interaction the 
width seems to increase explosively for large \po | . At even 
higher |po|, however, the opening of phase space is sup- 
pressed by the NJL cutoff and the width decreases again. 
Physically the most interesting area lies in the energy 
range < po < ^^^^ since that is the region of the pop- 
ulated quark states. All states above the Fermi energy 
as well as the anti-quark states at negative po are unoc- 
cupied (no holes in the Dirac sea). The structure of the 
spectral function is dominated by the on-shell peaks of 
the quarks and anti-quarks at po = |p| and pq = — |p|. 
The peaks get broader when the coupling is increased. 
Strength is distributed away from the peaks to the off-shell 
regions, the width of the peaks increases from 0.1 — 1 MeV 
to 10 MeV. The width seems to scale with the coupling 
strength (resp. G^, cf. eqs. (||)) while the general shape 
remains unchanged. 

Fig. ^ shows the width and the spectral function for 
the two chemical potentials ujp = 0.268 GeV and ujp = 
0.387 GeV using the regular coupling strength. The effect 
of the higher density is comparable to increasing the cou- 
pling. The width approximately scales with the chemical 
potential while the shape remains unchanged. This leads 
again to a broadening of the peaks of the spectral function. 

In the momentum distribution of nucleons in nuclear 
matter a depletion of the occupation probabilites by about 
10% is seen [l|. The resulting high energy tail is taken 
as a universal sign of short-range correlations. We show 
the momentum distribution of the quarks for the different 
coupling strengths in Fig. ^(a). At the lowest couphng a 
depletion of only 0.1% is found. For the coupling twice 



4 Summary and Outlook 

Based on a successful model for nuclear matter we have 
presented a simple but fully self-consistent approach to 
the spectral function of quarks in quark matter. It uses 
a pointlike interaction and goes beyond the quasiparticle 
approximation. The numerical results indicate that the 
influence of short-range correlations is rather small com- 
pared to nuclear matter. This finding might be an artifact 
of the present model, the NJL model with vacuum para- 
meters in the Born approximation. Furthermore we have 
not considered broken symmetries. However, the calcula- 
tions have shown that the model is technically feasible and 
suitable for further development. 

At present, the model is extended in two ways. First, 
a more realistic phase with broken symmetries - the chi- 
rally broken phase - is considered. In this phase large con- 
stituent quark masses must be taken into account and 
the analytic structure of the spectral function (||) is more 
complicated. Second, a new class of diagrams is incorpo- 
rated into the coUisional self-energies (|). This series of di- 
agrams is shown in Fig. |^. The lowest order contribution 
is again the Born diagram. The full series can be inter- 
preted as a dynamically generated meson that couples to 
the quarks ^M- First estimates show that this extension 
will increase the effective quark coupling and should lead 
to significantly higher widths of the spectral function. In 
addition it will be possible to investigate the properties of 
the self-consistently calculated pions and sigmas. 
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